FLOWING MAPS TO MINIMAL SURFACES: EXISTENCE AND 
UNIQUENESS OF SOLUTIONS 

MELANIE RUPFLIN 



Abstract. We study the new geometric flow that was introduced in [11] that evolves 
a pair of map and (domain) metric in such a way that it changes appropriate ini- 
tial data into branched minimal immersions. In the present paper we focus on the 
existence theory as well as the issue of uniqueness of solutions. We establish that a 
(weak) solution exists for as long as the metrics remain in a bounded region of mod- 
uli space, i.e. as long as the flow does not collapse a closed geodesic in the domain 
manifold to a point. Furthermore, we prove that this solution is unique in the class 
of all weak solutions with non-increasing energy. This work complements the paper 
[11] of Topping and the author where the flow was introduced and its asymptotic 
convergence to branched minimal immersions is discussed. 



1. Introduction 

Let M be a smooth closed orientable surface and let (N, Gn) be a (fixed) closed smooth 
Riemannian manifold of arbitrary dimension that we view as being isometrically im- 
mersed in M. K for some K £ N. 

For g a Riemannian metric on M and a map u : (M, g) — > (N, Gn) the Dirichlet energy 
is defined as 



E{u,g):= \ f \d 

1 J M 



u\ 2 dv g . 

M 



We remark that (u, g) is a critical point of E if and only if u is harmonic and weakly 
conformal, i.e. a branched minimal immersion or a constant map. In the present paper 
we establish the existence theory for the natural gradient flow of E (considered as a 
function of both the map and the domain metric) which was introduced in [11]. We refer 
to this joint paper of Topping and the author for the construction and the geometric 
background of this flow, but for convenience here recall the main points that led to the 
definition in [11]. 

We consider the negative gradient flow of E considered as a function of both the map 
and the domain metric, but taking into account the symmetries of E, that is the in- 
variance under conformal variations of the domain as well as under the pull-back by 
diffcomorphisms applied simultaneously to the metric and the map component. That is 
we consider E and its gradient flow on the set 

A = {[(u,g));g€M c ,u€C°°(M,N)} 

of equivalence classes where we identify (it, g) ~ (no/, /*<?) for smooth diffcomorphisms 
/ : M — > M homotopic to the identity. Here Ai c stands for the set of smooth metrics of 
constant (Gauss-)curvature c = 1, 0, —1 for surfaces of genus 7 = 0, 1 respectively 7 > 2, 
with unit area in case 7 = 1. 
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The tangent space of M. c splits orthogonally into a horizontal part consisting of the real 
parts of holomorphic quadratic differentials and a vertical part along the fibers of the 
action of diffcomorphisms on M c , i-e. the space of Lie-derivatives of the metric, compare 
Lemma 2.5 below. This canonical splitting allows us in [11] to represent solutions of the 
i 2 -negative gradient flow of E on A by the solutions of the system 

(1.1a) d t u = T g (u) 

(1-lb) J = ^Re(Pf($(u >ff ))). 

Here T g (u) = tr ff (V g (du)) = A g u + A g (u)(Vu, Vit), A the second fundamental form of 
TV <-} M. K , denotes the tension field of u : (M,g) — > (N,Gn) and $(u,g) stands for the 
Hopf-diffcrcntial, i.e. the quadratic differential given in conformal coordinates z — x + iy 
of (M,g) as g) — tpdz 2 for (j> — \u x \ 2 — \u y \ 2 — 2i(u x ,u y ). Furthermore denotes 
the L 2 -orthogonal projection from the space of quadratic differentials onto the finite 
dimensional subspacc of holomorphic quadratic differentials on (M,g). Finally r\ > is 
a free coupling constant related to the choice of L 2 -metric on A. 

As the main result of this paper we prove the following existence and uniqueness theorem 

Theorem 1.1. To any given initial data (uo,<7o) G C°°(M, TV) x A4 C there exists a weak 
solution (u,g) of (1.1) defined on a maximal interval [0,T) ; T < oo, that satisfies the 
following properties 

(i) The solution (u,g) is smooth away from at most finitely many singular times 
Ti G (0,T) at which 'harmonic spheres bubble off'. More precisely as t Tj 
energy concentrates at a finite number of points S(Ti) C M and suitable rescalings 
of the maps u(t) around points in S(Ti) converge as t Ti to (a bubble-tree of) 
non-trivial harmonic maps from R 2 U {oo} = S 2 to N. 
(ii) As t —¥ Ti the maps u(t) converge weakly in H 1 and smoothly away from the 
set S(Ti) to a limit u(Ti) £ H 1 (M, N). Furthermore, the metrics g(t) converge 
smoothly to an element g(T) G M c ; in fact, the flow of metrics is Lipschitz- 
continuous with respect to all C m metrics on A4 C across singular times. 
(Hi) The energy t i— > E(u(t), g(t)) is non-increasing. 

(iv) The solution exists as long as the metrics do not degenerate in moduli space; 
i.e. either T = oo or the length £(g(t)) of the shortest closed geodesic in {M, g(t)) 
converges to zero as t /T. 

Furthermore, the solution is uniquely determined by its initial data in the class of all 
weak solutions with non-increasing energy. 

Definition 1.2. We call (u,g) G H} oc (M x [0,T),2V) x C°([0,T), M-i) a weak solution 
of (1.1) if u solves (1.1a) in the sense of distributions and if g is piecewise C 1 (viewed 
as a map from [0,T) into the space of symmetric (0,2) tensors equipped with any C k 
metric, k G N) and satisfies (1.1b) away from times where it is not diffcrcntiablc. 

We remark that the assumption on the initial data in Theorem 1.1 can be weakened to 
Mo G H 1 (M, N) with the resulting solution being smooth away from finitely many times, 
possibly including T\ = 0. 

On intervals where the obtained solution (it, g) is smooth, the energy decays by 

(1.2) ±E(u(t),g(t)) = - fjr g (u)\ 2 dv g ~ ^\\Re[P g H (^(u, g))]\\l 2{MMt)) , 

so that if T = oo, both the tension field as well and the holomorphic part of the Hopf- 
diffcrcntial converge to zero as t — > oo suitably. In the joint paper [11] of Topping and the 
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author wc indeed prove that if the metric does not degenerate even as t — > oo, then the 
full Hopf-diffcrcntial (sub) converges to zero, resulting in a limit that is both harmonic 
and weakly conformal and thus, if non-constant, a minimal immersion away from at most 
finitely many branch-points [6]. More precisely, in [11], wc prove 

Theorem 1.3 ([11], Thm 1.4). In the setting of Theorem 1.1, if the length £(g(t)) of the 
shortest closed geodesic of (M,g(t)) is uniformly bounded below by a positive constant, 
then there exist a sequence of times U —> oo and a sequence of orientation-preserving 
diffeomorphisms fyiM—} M such that 

f*g(ti) -> g and u(U) o fa ->• u 

converge to a metric g S Ai c and a branched minimal immersion u or a constant map. 
Here the convergence of metrics is smooth, while the maps converge weakly in H 1 (M, N) 
and strongly in Wy'*!(M\S) for any p € [l,oo) away from a finite set of points where 
energy concentrates. 

For suitable initial data, such as incompressible maps, a degeneration of metrics can 
be excluded so that the flow (sub) converges (up to reparamctrisations) to a branched 
minimal immersion. In [11] we thus recover the well known results on the existence of 
branched minimal immersions with given action on the level of fundamental groups of 
Schoen-Yau [14] and Sacks-Uhlcnbcck [13] with a flow approach. 

Solutions of (1.1) that degenerate in moduli space will be analysed in a forthcoming 
paper [12] by Topping, Zhu and the author. 

Remark 1.4. For surfaces of genus less than two the structure of the flow (1.1) is 
simplified considerably and the existence of solutions is known; for spheres the space of 
holomorphic quadratic differentials is trivial so (1.1) reduces to the harmonic map flow 
of Eells and Sampson [3] for which existence of global weak solutions was proven in the 
seminal paper of Struwe [15]. For maps from a surface of genus 1 it is shown in [11] that 
(1.1) agrees with a flow that was introduced and studied by Ding, Li and Liu in [1]. In this 
special case the flow of metrics is reduced to two scalar ODEs for parameters describing 
a global horizontal submanifold of the space of metrics. Furthermore, the completeness 
of Teichmuller space prevents a degeneration of the metric at finite times, leading to the 
existence of global (weak) solutions for all initial data as obtained in [1]. 

In this paper we thus focus on the analysis of the flow from general surfaces of genus 
7 > 2. 

Outline of the paper 

The paper consists of three main parts. In the first section we study the properties of 
horizontal curves, i.e. curves that move in the direction of the real part of holomorphic 
quadratic differentials. Using ideas from Teichmuller theory, we obtain strong estimates 
for all horizontal curves, and thus in particular for the metric component of the flow, 
under the sole condition that we stay away from the boundary of moduli space. 

In the second section we prove the existence of solutions as claimed in Theorem 1.1. First 
we obtain short-time existence of smooth solutions based on the properties of horizontal 
curves derived in the first section. In a second step we then analyse the possible finite 
time singularities of the flow. On the one hand, we prove that the only way for the 
metric component to become singular is by a degeneration in moduli space. On the other 
hand, we obtain that as long as the metric component remains regular, the behaviour of 
solutions to (1.1a) is similar to the one of solutions of the harmonic map flow as described 
by Struwe in [15]; the singularity is caused by the bubbling off of harmonic spheres and 
the flow can be continued past the singular time by a weak solution. 
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Finally we consider the question of uniqueness. We show uniqueness not just for solutions 
of (1.1) satisfying properties (i)-(iii) of Theorem 1.1 but in the general class of weak 
solutions with non-increasing energy. This represents the analogue of the uniqueness 
results [4] and [5] of Freire for the harmonic map flow. 

Remark 1.5. For general curves within Ai-i, satisfying an L 2 bound on the velocity 
such as (1.2), singularities can form without the metrics degenerating in moduli space. 
For the flow 

(1.3) d t u = r g (u), g = ^Re($(« jfl )), 

which we would obtain if we were to consider the gradient flow of E without taking 
into account the symmetries, wc thus would not have a characterisation of the maximal 
existence time of solutions as statement (iv) of Theorem 1.1. For (1.3) we thus could 
not expect to obtain the global solutions needed to evolve pairs {u,g) to critical points 
of the energy, i.e. to branched minimal immersions, even for incompressible initial data. 

Acknowledgements: The author thanks Peter Topping for valuable discussions. This 
work was partially supported by The Leverhulme Trust. 

2. Horizontal curves 

We consider general horizontal curves, that is curves moving in the direction 

±g = Re(*(f)) 

of holomorphic quadratic differentials *f?(i) = ip(t)dz 2 on (M,g(t)), z = z{t) a complex 
coordinate on (M,g(t)). Key for the analysis of such curves is a good understanding of 
the dependence on the metric g £ M.-\ of the horizontal space 

H(g) := {Rc($) : $ = (j)dz 2 holomorphic quadratic differential on (M,g)} 

and of the corresponding L 2 -orthogonal projection Pg ■ What we essentially need is a 
quantified version of the idea that a smooth variation of the metric leads to a smooth 
variation of the complex structure, which in turn results in a smooth change of the space 
of holomorphic quadratic differentials and of P^ . 

We remark that there are several equivalent points of view that one can take to study 
horizontal tensors and curves as well as to study the flow (1.1). Here we follow the dif- 
ferential geometrical approach to Teichmiiller theory as presented in the book of Tromba 
[17]. We view the space of horizontal tensors as a subspace H(g) of the space Sym 2 (M) 
of all real symmetric (0, 2) tensors of class L 2 , with H(g) characterised by 

H(g) = {h e Sym 2 {M) : tt g (h) = and 5 g h = 0}, 

S g the divergence operator (induced by the Levi-Civita connection V g ). We then consider 
the projection 

P g : Sym 2 (M) -»• H(g) 
that is orthogonal with respect to the L 2 (M, <?)-inner product 

(k, h) L 2(M,g) '■= I g 13 g lm k a h jm dv g . 

J M 

This projection P g , for which we shall derive an explicit formula later on, is related to 
the projection P^ from the space of quadratic differentials to the space of holomorphic 
quadratic differentials by 

(2.1) Pg{Re^dz 2 )) = Re(P g H (^dz 2 )) 

for any quadratic differential ipdz 2 on (M,g). 
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We first remark that the set of hyperbolic metrics M.-\ with smooth coefficients (in the 
given coordinate charts) is not a manifold. On the other hand for any s > 3 the set 
MLi of hyperbolic metrics with coefficients in the Sobolev space H S (M) is a smooth 
submanifold of the (half-)space of all H s metrics on M, see [17] Theorem 1.6.1. 

We shall think of the projection P g as a map from this Banach manifold M. s _ l into the 
space L(Sym 2 (M) ,T Ai'Li) of linear functions mapping symmetric (0, 2)-tensors into the 
tangent bundle of Ai s _ 1 and prove that it is locally Lipschitz. 

Proposition 2.1. For any smooth hyperbolic metric go 6 M.-\ and every s > 3 
there exists a neighbourhood W of go in the Banach manifold M s _i and a constant 
C = C(go,s) < oo such that the following holds true: For every tensor k € Sym 2 (M) 
and every curve g £ C ,1 ([0, T), ML\) contained in W we have 

(2-2) \\P g (k)\\ H s<C-\\k\\ L2(Mig) 

and 

(2-3) ll^(t)(*)H*' <C-\\±g(t)\\H. ■ \\k\\ L *(M, g y 

Here and in the following the Sobolev norms || • \\h s are to be computed in fixed local 
coordinate charts of M. 

Based on this local statement about the projection P g , we then derive the following result 
for horizontal curves contained in compact regions of moduli space. 

Proposition 2.2. For every e > and every s > 3 there exists a number 9 = 8(e, s) > 
such that the following holds true. Let go £ M. s _i be any hyperbolic metric of class 
H s for which the length £(go) of the shortest closed geodesic in (M,go) is no less than 
e. Then there is a number C = C(go,s) < oo such that for any horizontal curve g € 
C 1 {[0,T),M ! L 1 ) with g(0) = g and of L 2 -length 




we have 

(2-4) \\^g(t)\\ H s<C\\^g(t)\\ L2(MMt)) forevewt&[0,n 

For tori the corresponding result is obtained as a consequence of the existence of a 
smooth global horizontal slice, i.e. of a finite dimensional smooth submanifold of A4o, 
parametrised over Teichmiillcr space, whose tangent space at each point is horizontal 
and which thus contains all horizontal curves passing through go- 

While for surfaces of genus 7 > 2 the space of horizontal tensors H(g) is still finite 
dimensional, dmis.(H(g)) = 67 — 6 by the Ricmann-Roch theorem, the distribution g i-> 
H (g) is no longer intcgrablc, compare [17], section 5.3, so Proposition 2.2 cannot be 
reduced to a statement about curves on a finite dimensional manifold. 

Proof of Proposition 2.1. We prove Proposition 2.1 in two steps; we show first that esti- 
mates of the form (2.2) and (2.3) hold true for metrics contained in a so called slice and 
then in a second step pull-back these estimates to give the claim of Proposition 2.1 for 
general metrics in a neighbourhood of go- To do so we make use of ideas from Teichmiillcr 
theory as explained in the book of Tromba [17], chapter 2. 

So let go e M~i be any given metric and let s > 3 be fixed. Following [17] we define a 
small slice around go by 

(2.5) S := {g = p(h) ■ (g + h) : h 6 U C H(g )} C M-i 
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for U = U (go, s) a suitably small neighbourhood of G H(go) chosen later on. Here the 
function p(h) : M — >• M. is to be chosen such that p(h) ■ (go + h) has constant curvature 
— 1 and is uniquely determined by this property according to Poincare's theorem. 

The key feature of this finite dimensional submanifold of Ai s _ 1 is that it provides a local 
model of M a _ x /Vl +1 , with V S Q +1 the set of £P +1 -diffeomorph isms that arc nomotopic to 
the identity 

Theorem 2.3 ([17], Thm 2.4.3). For any number s > 3, any g G M-i and S = S(go, s) 
a sufficiently small slice around go, there are neighbourhoods W C WL^ of go and V C 
2?q +1 of id for which the map 

SxV3(g,f)^ f*g G W 

is a diffeomorphism. 

For a proof of this theorem as well as for further insight into Teichmiiller theory we refer 
to the book of Tromba [17]. We remark that the above result remains valid if we replace 
the slice S by a smaller slice defined by (2.5), for appropriate new neighbourhoods of id 
in T>q +1 and of go in Ai s _i, but that the theorem does not give the existence of a uniform 
slice for which the statement is valid for all numbers s > 3. We furthermore stress that 
the theorem demands that the metric go is not only in M s _ 1 but smooth; this in turn 
implies that all metrics contained in a small slice S are smooth and thus satisfy stronger 
estimates than just H s bounds, in particular 

Lemma 2.4. For a sufficiently small slice S around go G M—i there exists a constant 
C = C(s,go) < oo such that for all metrics g\.2 G S 

(2-6) 1 1 51 - ffallff'+i < C- ds(gi,92) ■ 

Here we denote by ds the H s metric on S, i.e. consider S as a submanifold of the Banach 
manifold M. s _ v 

Apart from the finite dimensionality of H(go), and thus of S, the essential observation 
leading to the above estimate is that the conformal factor p(h) can be characterised as 
the unique solution of an elliptic PDE, compare [17] section 1.5, leading to a smooth 
dependence of p(h) on h G U. 

Based on these stronger estimates on elements of the slice, we can analyse the dependence 
of P g on g G S using an explicit formula for P g that we shall derive now. 

We first recall the following canonical splitting of the tangent space T g M s _ l into the 
horizontal and vertical space, sec Theorem 2.4.1 of [17]. 

Lemma 2.5. For any g G A4-i the tangent space T g M s _ 1 splits L 2 -orthogonally into 
H(g) and the space {Lxg} of Lie-derivatives. More precisely, given any k G T g M s _\ 
there is a unique vector field X (of class H s+1 ) such that 

tr g (k — Lxg) = and S g (k — Lxg) = 

and X can be characterised as the unique solution of the elliptic PDE 

(2-7) Sgd* g X = Sgk, 

6*X = —Lxg the L 2 (M, g)- adjoint of S g . 

In order to define the orthogonal projection of a general symmetric (0, 2) tensor k onto 
the horizontal space H(g), we first map k onto an element of T g A4 s _ 1 using 
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Lemma 2.6. For any g G M.-i and any symmetric (0,2) tensor k of class H s there 
exists a unique function fi G H S (M,S.) such that 

k-n-geT g M s _ v 
The function \i is characterised as the unique solution of the equation 

(2.8) -A gl i + 2 l x = 2DR{g){k), 

R{g) the Gauss curvature of(M,g). 

Given any g G M-i, we now claim that the orthogonal projection P g : Sym 2 (M) — > H(g) 
is given by 

(2-9) P g (k) :=k- v{k,g) ■ g - L X (k-n(k,g).g, g )9 

where X(-) and /z(-) stand for the corresponding solutions of (2.7) and (2.8). Indeed, 
Pg\n(g) = id and for general k G Sym 2 (M) the tensor given by (2.9) is well defined and 
divergence- as well as trace-free with respect to g, i.e. an element of H(g). Furthermore, 
k — P g (k) stands orthogonal to any h G H(g) as 

(h, k - P g {k)) L 2(M,g) = {h,(J,- g)h^ + (h, L x g)^ = fJ> • tr g (h) dv g + (h, -S*X) L 2 

Jm 

= -(Sgh,X) L 2 =0. 

To analyse the dependence of P g on g we now use that X and /i are characterised by 
elliptic PDEs for which the following uniform estimates apply 

Lemma 2.7. Let s > 3 and go G M-i be given and let S = S(go,s) be a sufficiently 
small slice. Then there exists a constant C = C(s,go) < oo such that the following 
claims hold true for every g G S. For every vector field Y there is a unique solution of 
the equation 

(2.10) 5 a 5* g X = Y 
and for any < I < s + 1 we have 

\\X\\ Hl <C-\\Y\\ Hl - 2 . 
Similarly, the unique solution fj. of 

(2.11) -A g ^ + 2iJ = feH l - 2 (M 1 R) 1 
satisfies 

<C-\\f\\ H i- 2 , 0<l<s. 

We remark that the occurring Sobolev norms with negative exponent are to be under- 
stood as the norms of the coefficients in the dual spaces H~ k (il) = (ifg(f2))*, il C M 2 . 

The reason why the solution X of (2.10) is unique is that we work on a surface that has 
negative curvature. Thus the kernel of S g 8*, which agrees with the space of Killing- fields, 
is trivial, see e.g. [7], Thm. 5.3. Elliptic regularity theory combined with the Fredholm 
alternative theorem then immediately gives the estimates for each individual g G S. 
These estimates arc indeed uniform since all metrics in S are contained in a small (H s ) 
neighbourhood of go. 

We can now give the proof of Proposition 2.1, first for metric contained in the slice. 

Let go G M-i, s > 3 and let S = S(go, s) be a small slice as defined above. Combining 
the elliptic estimates of Lemma 2.7 with (2.9) and the bounds on g given in Lemma 2.4, 
we find that for every < I < s and every k G Sym (M) 

\\Pg{k)\\w <\\k\\ w +c\M w +c\\x\\ Hl + 1 

< \\k\\ Hl + C(\\DR(g)(k)\\ Hl - 2 + \\6 g k\\ H i-i) < C\\k\\ Hl . 
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Here and in the following we crucially use that Lemma 2.4 gives bounds on s+1 derivatives 
of g so that we may estimate the H s and not just the H s ~ l norm of Lie-derivatives Lxg- 

Similarly, given any C 1 curve g in the slice, we differentiate the corresponding equations 
(2.7) and (2.8) characterising X and /i. This leads to elliptic PDEs of the form (2.10) 
and (2.11) for jfX(t) and ^/^(t). Applying Lemma 2.7 and making use of the bound 
II 37.9ll,ff s + 1 < C ' Ig^sll-H" 3 of Lemma 2.4, we obtain 

(2-13) \\^P gl t)(k)\\H><C\\±g\\ H .-\\k\\ Hl 

for any (sufficiently smooth) tensor k G Sym 2 (M) and any < I < s. 

In order to establish the estimates (2.2) and (2.3) claimed in Proposition 2.1 we now 
need to prove that the two estimates (2.12) and (2.13) obtained above remain valid with 
the H l norm on the right hand side replaced with the L 2 norm. We use 

Claim: There exists C < oo such that for all g G S and all h G H (g) 

\\h\\ H * < C\\h\\ LHM , g) . 

Proof of Claim: The estimate trivially holds true for g = go (or indeed for any one fixed 
metric) since H(go) is a finite dimensional space of smooth tensors. For general g G S 
we can parametrize H(g) over H(go) by restricting the projection P g onto H(go). Using 
estimate (2.12), we then get 

(2.14) \\P g {k)\\ H * < C\\k\\ Hs < C\\k\\ L 2 for every k G H(g ),g G S, 

with || • \\ L 2 denoting one of the equivalent L 2 (M,g) norms, g G S, say || • \\L 2 (M, go )- 

On the other hand, integrating (2.13) for I = along a suitable curve of metrics con- 
necting go to g and making use of the fact that P go \H(g a ) = we obtain that for any 
k G if (so) 

\\k\\ L 2 < \\P g {k) - k\\ L , + \\P g (k)\\ L 2 < Cd s (g,g ) ■ \\k\\ L , + \\P g (k)\\ L 2 

<ip|| i2 + ||P g (fc)|| i2 

provided the slice is chosen small enough. Combined with estimate (2.14) this implies the 
claim for tensors in the image P g (H(g )) C H(g) which must agree with H(g) because 
P g\H( go ) is injective and dim(iJ(g)) = dim(i? (g )) . 

Combining this claim with the estimate (2.12) for I — we have thus proved the first 
claim (2.2) of Proposition 2.1 for general tensors k G Sym 2 (M) and for metrics g G S in 
the slice. 

To obtain an improved version of (2.13), we write 

p 9 (t)(k) = P g (t)(P g (t )(k)) + P m (P g{t) (k) - P g{t0 ){k)) 

and estimate the derivative of the right hand side at t = to. Estimate (2.13), applied 
first for I = s and then for I = 0, combined with the estimate (2.2) we just proved then 
implies that for any k G Sym 2 {M) 

\\(j t P g (t)(k))(to)\\ H « < \\±g\\ H > ■ + \\P g (t o) (j t P g{ t)(k))\\H< 

< C\\±g\\ H . ■ \\k\\» + \\±P g (t){k)\\v> < c W^ t 9\\ H . ■ \\k\\ L *. 

This completes the proof of Proposition 2.1 for metrics g contained in the slice. We now 
pull back these estimates to the full H s neighbourhood W of go given by the slice-theorem 
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2.3. The key observation allowing us to do so is that the projection onto the horizontal 
space commutes with the pull-back 

f*P g (k) = P ra (f*k). 

Thus, given a C 1 curve g in W, we write it (uniquely) in the form g(t) = f(t)*g(t), for 
/(*) GFC £>o +1 and g(t) e S and recall that ||gj<7||ff* and \\-^f\\ H s+i are controlled 
by ||gj3||irs, see Theorem 2.3. Indeed, since the diffeomorphisms / are contained in 
a neighbourhood of the identity, also ||^/ _1 ||_ffs+i is bounded in this way. Applying 
estimates (2.2) and (2.3) for g e S, we thus find 

\\p B (k)\\B' = um((rTfc)ii^ < c ■ wPgdr'ym*' < c\\k\\„ 

as well as 

\\j t p g{t) (k)\\H» < c ■ \\^ t f\\ H » + i ■ wmr'rmiH* + c\\±(p m (mr x yk)\\H. 
< c ■ (ll^/ll^ + 11^511^ + ll^r W) • \\k\\ L * 

<C-\\±g\\ H ,-\\k\\ L * 

for any tensor k e Sym?(M) and any curve in W as claimed in Proposition 2.1. □ 



Proof of Proposition 2.2. For any number s > 3 we define a function 9 : M. s _ 1 —> [0, oo] 
as follows. For any metric go € A4 S _ 1 we let 0(go) be the supremum of all numbers 
6 > such that there exists a number C < oo for which estimate (2.4) holds true for 
all (piecewise) horizontal curves in A4 S _ 1 of length Lj l 2(g) < 9 and with g(0) = go. We 
stress that both this constant C, as well as the constant in Proposition 2.2, are allowed 
to depend on the metric go- 

We first claim that the function 9 is strictly positive for all smooth metrics. So let 
go e A4-i and let W be the neighbourhood of 170 m WL± for which Proposition 2.1 
applies. Writing the velocity of any horizontal curve as -^g = Pgi^g) and applying 
Proposition 2.1 we find that 



\^9\\ H s<C\\±9h* iM , g) 



for as long as the curve is contained in W . But W is an H 8 neighbourhood, so this 
estimate implies that any curve of small enough L 2 length and with g(0) = go is fully 
contained in W and thus that indeed 9(go) > 0. 

Secondly, we observe that 9 is invariant under the pull-back by diffeomorphisms. More 
precisely let T> s+1 be the set of all diffcomorphism of class H s+1 (not necessarily homo- 
topic to the identity). Then we claim that for any g e M s _ l and any / e V s+1 

d(fgo) = %o). 

Indeed, pulling-back any horizontal curve g in A4 S _ 1 by a fixed diffcomorphism / e T> a+1 
results in another horizontal curve of the same L 2 -length and with velocity bounded 
by ll^(/*5(i))||-ff s < C • ll^.9(0llff s j with C < 00 a constant depending on /. But we 
defined 9(g) asking only for an estimate of the form (2.4) to be satisfied for some constant 
C < 00, allowed to depend on the considered metric, so the claim follows. 

We conclude that 9 induces a positive map 9 on moduli space A4-i/T> and now want to 
prove that this function is continuous with respect to the Wcyl- Peterson metric dwp- 
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We recall that the length of a C 1 curve [g] in moduli space (with respect to the Weyl- 
Peterson metric) is given by 

Lwp{[g\) = 2 L L*(g) 

for g a 'horizontal lift' of [g], that is a horizontal curve g G M.-i with [g(t)] = [g(t)] for 
each t. 

Given any two points {g±\ and [32] in Ai^i/'D we now claim that 

%2])>%i])-2-d V y P ([ ffl ],[ ff2 ]), 

and thus switching the roles of [g\] and [52] that 9 is Lipschitz continuous on moduli space 
(M.-1/V, dwp)- So let 8 > be any fixed number and choose a (piecewise) horizontal 
path g of L 2 -length less than 2 • <iiyp([<?i], [92]) + 8/2 that connects a representative /*<?i 
of [gi] with g 2 . Let now g be any given (piecewise) horizontal curve with g(0) = gi and of 
length Lui{g) < 6([gi]) — 2dwp([gi], [32]) — 8. Precomposing it with g we obtain a curve 
G of length L L 2 (G) < 0([#i]) - (5/2 = 9(f* gi ) -8/2 and with starting point G(0) = 
By definition of 8(f*gi), the estimate (2.4) is satisfied for the extended curve G and thus 
in particular for g itself, with a constant C depending on f*gi and possibly (5 but not on 
g. We obtain the claim since 8 > can be chosen arbitrarily small. 

Given any number e > we now consider the subset K e of moduli space consisting of 
the equivalence classes of smooth metrics with shortest closed geodesic of length no less 
than s. This set K e is compact by the Mumford compactness theorem, see e.g. [17], p. 75. 
As a positive and continuous function, is thus bounded away from zero uniformly on 
K e which implies Proposition 2.2 for smooth metrics. 

For non-smooth metrics g G M. 8 _i \ M-i, we finally obtain the claim of Proposition 2.2 
using the invariance of 6 under H s+1 diffcomorphisms as well as 

Lemma 2.8. Given any g G A4 S _ 1 there exists a smooth metric g G M-i and a diffeo- 
morphism f of class H s+1 such that 

/* — 
9- 

For the sake of completeness we provide a proof of this fact in the appendix. 

For most arguments in the rest of the paper the estimates of Proposition 2.1 and 2.2, 
controlling the L 2 -orthogonal projection in terms of the L 2 norms of the involved tensors, 
would be sufficient, though would in some cases lead to slightly weaker regularity results. 
For the proof of uniqueness of weak solutions carried out in section 4 it is however crucial 
that we can extend P g continuously onto the space of tensors with finite L 1 norm 

Lemma 2.9. For any go G M.-\ and any s > 3 there exists a neighbourhood W of go 
in A4 S _ 1 such that the following holds true. The map P g is Lipschitz- continuous as a 
map from W to the space of linear maps from (Sym 2 (M), ||||x,i) to the tangent bundle 
TM. s _i, i.e. there exists a constant C = C(go,s) < 00 such that for all 51,52 G W and 
k G Sym 2 (M) 

(2.15) \\P gi m\H° <C||A;|| L i and \\P gi (k) - P g2 (k)\\ Hs < C ■ d M s_ i (g 1 ,g 2 ) ■ \\k\\ L u 

We remark that there is no need to specify with respect to which metric g G W the L 1 
norm is computed as all metrics in W are equivalent. 

We prove these refined estimates on P g using the following consequence of Proposition 
2.1 
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Lemma 2.10. For any go 6 M.-1 and any s > 3 there exists a neighbourhood W of 
go in Mlj and a constant C < oo so that we can assign to each metric g in W an 
L 2 (M, g)-orthonormal basis {©^(ff)}^^ 6 °f H(.9) satisfying 

\\®°(9i) ~ ® j {92)\\H° < C ■ d Mti (g u g 2 ), g ij2 eW, j = 1 . . . 67 - 6 = dim(H(g)) 

Lemma 2.9 then immediately follows from P g (k) = ^2j{k, © 3 (g))L 2 (M.g)0"' (.<?)■ 

Proof of Lemma 2.10. Let go € A4 S _ 1 and let W be the neighbourhood of go given by 
Proposition 2.1. We fix any L 2 (M, po)-orthonormal basis 6 J (<7o), j = I...67 — 6, of 
H{go) and define 

e ( g ) : /yo<;,/„n. 

According to Proposition 2.1 this auxiliary family of tensors depends continuously on g, 
\H(9i) ~ ®i(92)\\ H . < C ■ d M '_M>9*) 

so that {0o} is a basis of H(g) provided the neighbourhood W is chosen sufficiently 
small. Furthermore, as the map asigning to each metric g the inner products 

g^ (el(g),e k Q (g)) L 2 {M!g} 

is also Lipschitz-continuous on W, so are the coefficients a\ of the orthonormal basis 
^(9) = Ei=i a i (d)® 0(9) °^ H(g) obtained by Gram-Schmidt orthogonalisation and 
thus the basis itself. □ 



3. Existence of solutions 



In this section we establish the existence of weak solutions to (1.1) satisfying the prop- 
erties claimed in Theorem 1.1, in particular existing for all times unless the metric com- 
ponent degenerates in moduli space. As a first step, we prove the following short-time 
existence result 

Lemma 3.1. For any initial metric go G M.-i and any initial map uo €E C°°{M,N) 
there exists a smooth solution (u,g) of equation (1.1) to initial data (u(0), g(0)) = (ito, <?o) 
defined on an interval [0, T), T = T(uo,go) > 0. 

Proof of Lemma 3.1. We first recall that the metric evolves by 

(3.1) J = ^Re(P g H = £p g (k(«,9)), 

where k(u,g) — Re($>(u, g)), compare (1.1b) and (2.1). 

To simplify notations and without loss of generality, we shall from now on consider the 
flow with coupling constant r\ = 2. We also remark that computing the variation 

^-E(u,g + sl)\ s= o = -](Re($(u,g)),l) L 2 for all I e Sym 2 (M) 
as 4 

in local coordinate charts, allows us to write the real part of the Hopf-differcntial in 

general (not necessarily conformal) coordinate charts as 

k(u,g) = Re($(u,g)) = 2u*G N - 2e{u,g)g, 

e(u,g) = \ |Vu| g = \g tj d Xl u ■ d Xj u the energy density. 

Using the results of the previous section we can consider equation (1.1) as a system 
consisting of a semilinear parabolic PDE coupled with a differential equation on a Banach 
manifold Ai 3 _ 1 that is defined by a locally Lipschitz continuous vector field. In such a 
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setting we obtain the existence of a classical solution on a short time interval using 
a standard iteration argument, which, for the sake of completeness, we outline in the 
appendix. Given any («o,<7o) G C 2 ' a (M, N) x M-i and any number s > 3 we obtain a 
solution 

(u,g) g C 2a ' Q ([0,T s ) x M,iV) x C 1 ([0,T),.Mi 1 ) 

of (1.1), defined on a maximal interval [0,T S ). This interval might a priori depend not 
only on (mo, <?o) but also on the Banach manifold M s _ 1 on which we solve (1.1b). Indeed, 
the key step needed to prove that the obtained solution (it, g) is actually smooth is to 
show that this is not the case. So suppose that for some 3 < si < S2 wc have T Sl ^ T S2 . 
Since classical solutions of (1.1) arc uniquely determined by their initial data, compare 
section 4, we remark that the two solutions obtained for the different values of s agree 
for as long as they both exist, that is until time T S2 < T Sl . Since the metric component 
is continuous (as a map into M. 8 2\) up to time T Sl there exists a number s > such 
that the length £(g(t)) of the shortest closed geodesic of (M,g(t)) is no less than e on 
the smaller interval [0,T S2 ]. Using the H s estimates of Proposition 2.2 this allows us to 
conclude that g is C 1 as a curve into M S J ± on the closed interval [0,T S2 ], compare with 
the proof of Lemma 3.2 below. 

Using Lemma 2.8, we then write g(T S2 ) £ M."2\ in the form g(T e2 ) = f*g(T S2 ) for an 
H S2+1 diffeomorphism / and a smooth metric g 6 M-\. Restarting the flow with the 
pulled-back initial data (u(T S2 ),g(T 32 )) = (u(T S2 ) o /-*, (/ -1 )*3(T S2 )) g C* 2 < a x M-i 
we obtain a solution {u,g) of (1.1) in C 2 ' 1,Q (M x I) x C 1 (/,A^'! 2 1 ) on a time interval 
/ = [T S2 ,T S2 +5). But equation (1.1) is invariant under the pull-back by diffcomorphisms 
applied simultaneously to both the map and the metric component and solutions of (1.1) 
are unique. Thus the pull-back of (u, g) by / is nothing else than our original solution 
(u, g) so that g is in C 1 ([0,T S2 + S),^^), leading to a contradiction. 

At this point we are now in a position to argue by a standard bootstrapping argument, 
using parabolic regularity theory to improve the regularity of it, as well as the explicit 
formula for P g given in (2.9) to analyse higher order time derivatives of g. We obtain 
that (it, g) is indeed smooth. □ 

We remark that the results of section 2 allow us not only to establish short-time existence 
of solutions to (1.1) but already give the following characterisation of the behaviour of 
the metric component at a singular time 

Lemma 3.2. Let (u,g) be a smooth solution of (1.1) defined (and smooth) on a maximal 
interval [0,Ti). Then one of the following three statements holds 

(i) T\ — co, or 

(ii) T\ < oo but as t /* T\ the metrics g{t) converge smoothly to a limit g{T\) G M.-\; 
indeed g can be extended to a Lipschitz continuous curve from the closed interval 
[0, Ti] into each of the Banach manifold Ai'Li, s > 3, or 
(Hi) the metrics degenerate in moduli space at a finite time Ti, i.e. lim^Xi ^(ffW) = 
0. 

Proof of Lemma 3.2. Assume that T\ < oo and that the length of the shortest closed 
geodesies in (M,g(t)) does not converge to zero 

limsup£(g(i)) > e > 0. 

Then given any number s > 3 we let 9 = 9(s, e) > be the constant of Proposition 2.2. 
We recall that according to the energy identity (1.2) the L 2 -length of the curve g is finite 
on intervals of finite length. We may thus choose to < T\ with £(g(to)) > e and close 
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enough to T\ such that L L 2(g\[ toTl - ) ) < 9. Proposition 2.2 then implies that g{t) is a 
Cauchy sequence in M s _ 1 and thus converges to a limit g(T\) in M s _ 1 as t /■ T\. Indeed, 
combining Proposition 2.2 with the energy identity (1.2) gives C 1 / 2 -H61der estimates in 
time for g considered as map into M. s _ x . Moreover, thanks to the uniform bound on the 
energy of u and thus on the L 1 norm of the Hopf-differential 

\\k(u,g)\\ L i<C-\\Vu\\ 2 L2 <C-E(u,g)<C-E(u ,g ), 

the improved estimates on P g stated in Lemma 2.9 give uniform bounds on || -^g{t) \\h s ■ 
Thus g is not only C 1 ' 2 but indeed Lipschitz continuous with respect to each H s metric 
on the closed interval [0, Xf]. □ 



We remark that the possibility of solutions degenerating in moduli space will be addressed 
in future work and that here we focus on the analysis of singularities of the second type, 
essentially due to the map component becoming singular. 

So let (u,g) be a smooth solution of (1.1) on a maximal interval [0,T). Assume that the 
metrics do not degenerate in moduli space as we approach the singular time and thus 
that g(t) —> g{Ti) £ M-x smoothly as t /*• T\. We remark that the evolution of the 
metric component is uniformly controlled, 

(3.2) < C\\k(u,g)\\ Ll < C ■ \\Vu\\ 2 L , < C ■ E 

for times in an interval of length 6 = d(g(Ti), s) > not just for the one solution (u, g) 
of (1.1) that becomes singular, but also for all solutions evolving from an initial data 
(u,g(Ti)) with energy bounded by Eq. Thanks to this strong bound on the metric 
component we can carry out the analysis of the map component of solution to (1.1) near 
singular times using methods familiar from the work of Struwe [15] on the harmonic map 
flow. Since our analysis closely follows the ideas of [15] we shall omit some details and 
calculations in the following presentation. We also remark that a similar argument was 
briefly outlined in [1] in the special case of maps from a torus. 

Notation: We let ji £ A'f-i be a fixed metric that should be thought of as a limiting 
metric of a solution of (1.1) at a singular time. Then unless indicated otherwise all 
occurring objects such as norms, operators (like A), integrals, balls and so on are to 
be understood as the corresponding objects on the fixed Riemannian surface (M,gi). 
Furthermore, we denote generic constants (allowed to change from line to line) by C 
in case they depend only on g\ and Eq and will indicate any dependence on additional 
quantities accordingly. 

Based on (3.2) we henceforth restrict our attention to solutions of (1.1) satisfying 

(3.3) ||.9i -g(t)\\ H . < ei 

for some fixed number s > 3 and a small e\ = £1(51, s) > 0, chosen in particular such 
that Lemma 2.9 applies on this A4 S _ 1 neighbourhood of g\. 

We first remark that the evolution of the local energy is controlled by 

Lemma 3.3. For solutions (u,g) of (1.1) satisfying (3.3) the following local energy 
bounds hold true for any point x S M and any radius < r < Ti n j 

E(u(t), B r/2 {x)) < 2E{u{0),B r (x)) + C-^ 

and 

E(u(t),B r (x)) > )-E{u{Q),B r/2 {x))-A [ f ip 2 \d t u\ 2 dvdt ~ c\. 

1 Jo Jm r 
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Sketch of proof. Given x £ (M,gi) and < r < r in j(M,gi) we let <p £ Cg° {B r (x) 7 [0, 1]) 
be a standard cut-off function, i.e. such that tp = 1 on B r / 2 (x) and \Vp>\ < ^. A short 
calculation shows that for a solution (u,g) of (1.1) 



= y ip 2 \dtu\ 2 dv — J p 2 dtu ■ A g ( t }udv 
Id /* 

^ 2 |9 t u| 2 dv+-- p 2 \Vu\ 2 g{t) dv + R(u(t),g(t)) 



(3.4) 

= / to 2 \r).n.\ 2 d.v 4- 

2dt . 

with an error term that is bounded by 

\R(u,g)\<(°+C\\±g\\c°)-E(u,B r (x)) + l f p 2 \d t u\ 2 dv. 
t dt 8 J M 

Since || gj<?||c° i s uniformly bounded, this estimate integrates to give an upper and a lower 
bound on J ip 2 \Vu\ g / t -. — f p 2 |Vu(0)| ff ^ ^ . Combined with the fact that < g < 2g for 
all g,g satisfying (3.3), we obtain the claims of Lemma 3.3. □ 



An important consequence of the previous calculation is 

Corollary 3.4. Suppose (u,g) is a smooth solution of (1.1) defined on a maximal interval 
[0, Ti) for which (3.3) is satisfied. Then for any Eq > the set of points 

S :={x £ M : ]im sup E(u(t), B R (x)) > e for all R > 0}, 

is finite. 



In fact #5 < E /sq, since energy concentrates near points of S not just along a suitable 
sequence tj /• T\ but indeed for all sequences t /• T\, compare (3.4). 

Away from the finite set S we control the map component of the flow using the following 
lemma which should be seen as the analogue of Lemmas 3.10 and 3.10' of [15] 

Lemma 3.5. There exists a number > depending only on g\ and Eq such that 
the following statement holds true. Let (u,g) be a smooth solution of (1.1) on an open 
interval (0, T) and assume that (3.3) is satisfied. Let M' C M be an open set such that 
there exists a number R > with 

(3.5) E{u{t), B R (x)) < e for all (x, t) £ M' x (0, T). 

Then the parabolic Holder-norms of u and its spatial derivatives (upto order s — 2) are 
bounded uniformly on the sets [t,T] x M' , t > 0, with bounds depending only on r, R, 
gi, T , s, M' and the energy bound Eq. 

Remark 3.6. If the initial map u(0) is smooth on a neighbourhood of M' then the 
above result can be extended to give bounds on the Holder norms of u\m' and its spatial 
derivatives on M' up to time t = 0, now with bounds depending additionally on it(0), 
compare with remark 3.11 and 3.11' of [15]. 

Remark 3.7. Because of the non-local nature of the projection operator P g these esti- 
mates on u\m> allow us to improve the regularity of g\m' from the a priori known C ' 1 
dependence on time only in case M = M'. For M' / Mwe can improve the bounds 
of Lemma 3.5 to give C l a bounds in time on u\m' and its spatial derivatives while for 
M = M', i.e. away from singluar times, a bootstrapping argument gives estimates on 
any C k norm (in space and time) of (u, g) in terms of the quantities specified in Lemma 
3.5. 



Proof of Lemma 3.5. For the proof of this lemma we follow largely the ideas of [15]. We 
make use of the well known interpolation estimate, see e.g. [2] 
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Lemma 3.8. There are numbers £q > and C < oo (depending on (M,gi) and the 
target manifold) such that for all maps u G H 2 {M,N), a bound on the local energy of 

E{u,B r {x)) <e 

implies an H 2 -bound of the form 

(3.6) / if 2 \V 2 u\ 2 dv < -zE(u,B r (x)) + C f p 2 |t(w)| 2 dv, 



as well as an estimate of 



(3.7) / v? 2 | Vu| 4 dv < CE{u, B r {x)) 



^E(u,B r (x)) + I ip 2 \r(u)\ 2 dv 



Here ip G Cq° (B r (x)) denotes a cut-off function. 

Let now (u,g), M' and R > be as in Lemma 3.5, let x G M' and choose a cut-off 
function tp G C^(Bft/ 2 (x)). We first remark that for e\ > sufficiently small, the 
pointwise bound \t 9 (u) — t(u)| < Ce±(\ V 2 u\ + Vu| 2 ) implies that (3.6) and (3.7) remain 
valid with r(u(t)) replaced by T g ^(u(t)) = dtu{t). 

As in [15] we now differentiate equation (1.1a) in time and multiply with (p 2 dtu. Af- 
ter carefully analysing all occurring terms, in particular the terms due to the time- 
dependence of the metric, we find 

~ f if 2 \d t u\ 2 + ftp 2 \Vd t u\ 2 < (1 + Cex) ■ f if 2 \Vd t u\ 2 + C fcp 2 \d t u\ 2 \\7u\ 2 



C(R)(l + \\^ \\ 2 C 2 + J \d, 



it*! 8 ' 



Since we know that the H s norm, and thus also the C 2 norm, of ^.g is uniformly bounded 
by a multiple of the energy, we obtain that for e\ = £i(</i) > chosen small enough 

(3.8) A J <f \ dtU \ 2 + J<p 2 \V8 t u\ 2 < C(R)(l + ^) + cJ<p 2 \d t u\ 2 |Vu| 2 . 
Using the Sobolev embedding W L 2 as well as Lemma 3.8, we find 

.2 in „.|2 ,„ ,2 ^ , a ,2 n _ / f 2 |_|4a/2 



C J ^ 2 \d t u\ 2 \Vu\ z <cy\d tU \ 2 \\ L 2-( j ^ 2 |Vu| 4 ) J 

<C(\\\7^\d t u\ 2 )\\Li + \\\dtu\ 2 \\ L i)E 1 /2 [C(R)+ J V 2 \d t u\ 2 } 1/2 
<J<p 2 |V9 tU | 2 + Ce ^- J V 2 \d t u\ 2 + C(R) ■ (^ + 1). 
Integrating the resulting estimate (3.8) over any interval [^1,^2] C (0,T) we thus obtain 
/ ip 2 \d t u\ 2 dv * 2 < Ce • sup / ip 2 \d t u\ 2 dv + C{R,T). 

After possibly reducing e = £o(Eq, gi) > so that the factor Ceo < |, we conclude that 
for any r > 

sup / \d t u\ 2 dv<2 inf f ip 2 \d t u(t)\ 2 dv + C{R,T) 

te[r,T) J B R/i (x) t&[0,r)J 



< 2— + C(R,T). 

T 
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Repeating the above argument for a finite cover of balls .B.R/4(a:i) of M' we obtain a 
uniform estimate of 

f \d t u(t)\ 2 dv<C for all t £ [r,T) 
Ju 

on a small neighbourhood U of M'. According to Lemma 3.8 this implies a bound on 
Jjj, |V 2 u(i)| dv on a slightly smaller neighbourhood of M'. Applying Sobolcv's embed- 
ding theorem we then conclude that for any exponent p < oo 

f \Vu{t)\ p dv < C p te[r,T). 

JU' 

Then as in [15] we think of (1.1a) as an inhomogcneous heat equation 

d t u- A g u = A g (u)(Vu,Vu) G L p (U' x [r,T)) 

allowing us to apply standard regularity results for parabolic equations, see e.g. [8], chap- 
ter VII; we get bounds in the parabolic Sobolev-spaces W 2 ' 1 , and thus in the parabolic 
Holder spaces C a , on sets M' x [t', T],t' > t. We finally obtain estimates on the Holder 
norms of spatial derivatives of u (up to order s — 2) by a standard bootstrapping ar- 
gument which relies on the strong bounds on the velocity of horizontal curves given in 
Lemma 2.9. □ 



Let now (u,g) be a smooth solution of (1.1) on [0,Ti) whose metric component does 
not degenerate and thus smoothly converges g(t) — > g{T\) =: gi G M-i as t /• T\. We 
first remark that the uniform bounds on the energies E(u(t)) < 2E(u(t),g(t)) < 2Eq 
combined with the fact that dtu G L 2 ([0,Ti) x M) imply that the maps u(t) converges 
weakly in H 1 (M) to a limit u(T{) as t T\. Additionally, Lemma 3.5 gives uniform 
Holder bounds on u and its spatial derivatives away from the finite set S of concentration 
points so that u(t) converges also in Cf£ c {M \ S). 

We now remark that any concentration of energy must be due to the so-called bubbling 
off of (at least) one harmonic sphere. Indeed, the analysis carried out in [15] (p. 578/9) 
remains unchanged as long as the local energy estimates and H 2 bounds used in [15] are 
replaced by Lemmas 3.3 and 3.8. We obtain the following: For any point x <G S there 
are sequences of times ti /* T±, radii r.; — > and points Xi — > xq with energies on balls 
around Xi of 

E(u(ti),B 2n (xi)) < e and E(u(U), B r .(xi)) > ce , c = c{M,g\)>Q 
and with tension satisfying 




The rescaled maps Ui(x) = u(exp x (rix),ti), defined on larger and larger subsets of R 2 
are then bounded uniformly in H 2 and subconverge (weakly in H 2 , strongly in W^' p , 
p < oo) to a non-constant harmonic map of finite energy that is defined on R 2 U{oo} = S* 2 , 
called a harmonic sphere or bubble. The amount of energy that concentrates near Xo, 
and that is consequently lost as we pass to the limit t Ti, is no less than e* = e*(N), 
the minimal energy of such a non-constant harmonic map from S 2 to the target. 

Finally, as in [15], we (weakly) continue the flow past any such singular time by restarting 
from initial data <?(Ti) G A4_i and u(Ti) £ H X {M,N) n C^ C {M \ S) as follows: let 
Ujfi G C°°(M,N) be a sequence of maps that converge to u(T{) in H 1 (M) as well as 
in Cfo c (M \ S) and let (uj,gj) be the smooth solution of (1.1) corresponding to initial 
data (iijfi, g(T\)) that exists at least on some interval [Ti,Ti + Sj) according to Lemma 
3.1. We remark that the metrics gj are uniformly Lipschitz continuous and thus that the 
estimates derived above can be applied on [Ti, Ti+min(5j-, 6q)) for a number So(g(Tx)) > 
independent of j. 



FLOWING TO MINIMAL SURFACES: EXISTENCE AND UNIQUENESS 



17 



We now choose r > such that sup x&M E(uj t o, B r {x)) < £o/4, which is possible due to 
the strong H ^convergence of the initial maps. Then the local energy estimates of Lemma 
3.3 imply that there is no concentration of energy and thus in particular no blow-up for 
any of the maps Uj, on a uniform interval / = [Ti,Ti + cr 2 ], c = c(gi) > 0, in the sense 
that 

E(uj(t), B r/2 (x)) < e a for all j € N, x G M, t e I. 

According to Lemma 3.5 as well as remarks 3.6 and 3.7 we thus obtain uniform C ' a 
estimates in time for the maps Uj (and their spatial derivatives) in every compact subset 
of Mx [Ti, Ti+cr 2 ]\(S x {Ti}). Away from the singular time, we furthermore get uniform 
bounds on all C norms of (uj , g.j ) in space-time. We conclude that a subsequence of 
(uj,gj) converges smoothly on M x (Ti,Ti + cr 2 ] to a pair (u,g) which solves (1.1) 
classically on (Ti , T\ + cr 2 ] and weakly on [Ti , T\ + cr 2 ] . This solution achieves the initial 
data (u(Ti), g(Ti)) in the sense that for t\Ti the maps u(t) converges to u(Ti) weakly 
in JT^M) and smoothly away from the set S while the metric component g is Lipschitz- 
continuous across the singular time. Since the energy of the approximating solutions 
(uj,gj) is no more than Eiiij^, g{T\)) — > E(u(T\), g(Ti)) 7 the extended weak solution 
(u,g) has non-increasing energy also across the singular time T\. In particular the total 
number of all singular points UiS(Ti) x {T;} of such a solution is bounded by MhM^MI _ 
After possibly repeating the above argument to analyse any further singularities, we thus 
obtain a weak solution satisfying the properties (i)-(iii) of Theorem 1.1 and existing for 
as long as the metrics do not degenerate in moduli space. 

4. Uniqueness of weak solutions 

We finally discuss the issue of uniqueness of weak solutions. We prove that the solution 
(u,g) of (1.1) constructed in the previous section is uniquely determined by its initial 
data, not only among all solutions satisfying the properties of Theorem 1.1, but in the 
natural class of all weak solutions with non-increasing energy. We remark that a further 
argument as carried out in [10] actually gives uniqueness under the weaker assumption 
that the total energy does not instantaneously increase by more than a certain quantum 
at any time. We also remark that it is necessary to impose restrictions on the evolution 
of the total energy in view of the possibility of reverse bubbling, see [16]. 

So let (itj, <?i)i=i,2 be two weak solutions of (1.1) defined on an interval [0,T) that evolve 
from the same intial data 

(ui,ffi)(0) = (« ,5d) = (w 2 ,3 2 )(0) G H\M,N) x M-i 

and assume that the total energies 1 1-> E(ui^(t), gi,2{t)) are non-increasing. Since 

I:={tE [0,T) : (u u9l ) = (u 2 ,g 2 ) on [0,t]} 

is trivially closed in [0,T), we need to prove that / is also open. 

Given any to £ I we recall that gi(t) — > <7j(io) in each M. s _ 1 and thus certainly uniformly 
as t \ to. Combined with the fact that Ui(t) — > u.i(to) strongly in L 2 and weakly in H 1 , 
we thus obtain 

E(ui(t ),gi(to)) > lim E( Ui (t), 9i (t)) = Hm E( Ui (t), gi (to)) > E( Ui (to),gi(t )), 

where we used the assumption on the evolution of the energy in the first step. Thus 
Ui(t) — > Ui(to) indeed strongly in H (M, go), which implies in particular that local ener- 
gies, say on balls, converge as t \ to- Choosing a finite cover of balls B r {xi), i = 1 . . . K, 
of (M,gi(to)) such that E(m(to), gi(to), B2 r {xi)) < en/2, we may thus choose S > so 
small that 

E(u l (t),g l (t ),B 2r (x l )) < e for t e [t ,t + S], i = 1,2. 
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Here we let Eq > be the constant of Lemma 3.8. 

It is now crucial to remark that on almost every time slice the functions Ui(t) weakly 
solves an almost harmonic map equation, that is an equation of the form r g v — f for a 
function / G L 2 and a metric g, here of course g = gi(t) and / = dtUi(t). Since any 
weak solution of such an elliptic equation is contained in the Sobolev space H 2 , see e.g. 
[10], Proposition 2.1, we may apply Lemma 3.8 on almost every time slice resulting in 
an estimate of 



for t G [to, to + 5) and i — 1,2. We can thus reduce the uniqueness statement in the 
general class of weak solutions with non-increasing energy to the following lemma whose 
analogue for the harmonic map flow was proven in [15] 

Lemma 4.1. Let (ui,gi) and (u 2 ,g 2 ) be weak solutions of (1.1) to the same initial data 
(iti, <7i)(0) = (u 2 ,g 2 )(0) and suppose that 



Then (iti,3i) = (u 2 ,g 2 ). 

Proof of Lemma 4-1- Using an open-closed argument as above it is enough to prove that 
the solutions agree on a possibly smaller interval [0, <5) , which we can chose in particular 
such that the metrics g\ t2 are contained in an H s neighbourhood of go = <7i(0) = 52(0) 
for which Lemma 2.9 applies. Here s can be chosen to be any fixed number s > 3. 

Notation: For the following computations we denote by || • \\lv the L p (M,g ) norm and 
by d(-, •) the metric on Ai'Lx respectively by [| • || the H s norm on TA4 S _ 1 . Furthermore, 
we use the short-hand notation of |VV| := max{|Vui| , VM2I} which, by assumption, 
is a function in L A {M x [0,T]) with L 2 norm on time-slices bounded by the energy, 
||W(t)||2 2(M) < C ■ (E(u 1 ,g 1 )+E(u 2 ,g 2 )) < CE(u ,g ). 

Subtracting the equations (1.1a) for the map components Ui we obtain that the difference 
iv = u\ — u 2 satisfies 

(4.2) d t w -A gi w= (A gi - A 92 )(u 2 ) + A gi (ui)(Vui, Vu x ) - A g2 (u 2 )(Vu 2 ,\7u 2 ) 

where A denotes the second fundamental form of the target N <^-t M . , A g (u)(Vu, Vm) := 
g % i A(u)(diU, dju). 

Following [15] we multiply equation (4.2) with w, integrate over the fixed surface (M,go) 
and estimate the resulting terms using Holder's inequality. This leads to 




(4.1) 



V Ul G L i {M x [0,T)), and V 2 Ui G L 2 {M x [0,T)) i = 1,2. 



Mh + \\Vw\\ 2 L2 <C-d( gi ,g )-\\Vw\\l 



(4.3) 



+ C ■ d( 9l ,g 2 ) ■ (\\Vw\\ L 2 + \\w\\ L 2 ■ (1 + ||W||i«)) 
+ C||V W || L2 • ||VV|| L 4 • \\w\\ L * + C||V^|||4 • ||«;|||4 

<{\ + C-t)\\Vw\\l 2 +C-d{ gi ,g 2 ) 2 

+ c(i + ||vF|||.)-|klli- 
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forV(t) = (||Vy(i)||l 4 + l)eL 1 ([0,T]). 



r(.9i - .92) = P gi (fc(wi,5i)) - P g2 (k(u 2 ,g 2 )), 



In order to estimate the distance of the metric components g\ and g 2 in terms of w wc 
recall that the evolution of the tensor g\ — g 2 is given by 

d 

k(u,g) = 2u*Gn — 2e(u,g)g. We have a pointwisc estimate of the difference of the 
involved tensors of 

|fc(ui, 5 i) - k(u 2 ,g 2 )\ < C ■ d(gi,g 2 ) • |W| 2 + C ■ \w\ ■ |W| 2 + C ■ \Vw\ ■ |W| . 

Remark, that any L 2 estimate of this tensor would involve integrals of the form J | X/V\ A |u>| 2 
and J \W\ 2 ■ \\7w\ 2 which are not controlled by the quantities of the left hand side of 
(4.3). It thus crucial at this point that the improved bounds on P g given in Lemma 2.9 
only ask for L 1 bounds on the involved tensors, allowing us to estimate 

^ 44 ^ ^d(gx,g 2 ) < C-d(g 1 ,g 2 ) ■ \\k(ux,gx)\\^ + C ■ \\k(ux,gx) - k{u 2 ,g 2 )\\ L i 

< C ■ d( 9l ,g 2 ) + C ■ \\Vw\\ L 2 + C ■ m 1/2 \H\ L i- 
GronvalPs lemma thus leads to an estimate of 



<%i,g 2 )(t) 2 <C- ( [ \\Vw(s)\\ LHM) ds) +c( f ^(s) 1/2 -|k( S )|U 2( A/) 

(4 5) 



2 

ds 



<t- I I \Vw\ 2 + C i)(s)ds- / w 2 , 

JO J M JO Jo J M 

which we insert into (4.3). Integrating the resulting estimate over time, we find 

f [ |Vw| 2 <C- f ^(a)da- sup \\w(s)\\ 2 L2 + Ct 2 f f \Vw\ 2 . 
Jo J Jo se[o.t] Jo J 

Since tp is integrable, we conclude that for all t sufficiently small, say t £ (0,£o); 

/ \w(t)\ 2 < - sup f \w(s)\ 2 . 
J z se[o,t] J 

Thus w must vanish identically on (0, to) so u\ = u 2 and g\ = g 2 as desired. □ 



A. Appendix 



A.I. Solving the equation on a fixed Banach manifold. 

Let (uo,go) £ C 2,a (M) x M-i, a > be given and let s > 3 be a fixed number. 
Here we outline an iteration argument that can be used to obtain a solution (it, g) G 
C 2 ^ a ([0,S) x M) x C 1 ([0,<5),.M1 1 ) of (1.1) for such initial data. 

For <5o = 5o(uo;50; s) > to be determined later, we extend uq to a constant in time map 
defined on M x [0, <5q) and define iteratively for i = 1 . . . 

• gi £ C^QO, M s _ x ) as the solution of ^.g 4 = P gi (k(ui-x, gt)) with ^(0) = g ; 

• Ui e C f2 ' 1 ' a (M,5 t ) as the solution of d t u.i = r gi (ui), Uj(0) = u , defined and 
smooth on a maximal domain M$ t := M x [0, Sf), 6i < Si-x- 

Here wc use the Lipschitz-continuity of the map P g on the Banach manifold Ai s _ 1 in the 
first step. Wc also remark that the equation for m is a semilinear parabolic equation so 
standard methods, see e.g. [9] Theorem 5.2.1, lead to the existence of a solution m of 
the above equation, defined on all of [0, Sf^) unless there is a blow-up in the gradient at 
some time Sf, < <5, < Si-x- 
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We claim that for So initially chosen small enough, the iterates are all defined on [0,<5q) 
and satisfy 

(A ^ hi+l - Ui\\c».i.-(M t0 ) ^ \\\ u * - u i-l\\*C^{M So ) 

- 9iVci([0,g*],M*-i) - c - s o\\ui + i ~ Ui\\* c2 ,i, a ^ Msa) , 

thus converging to a classical solution (u,g) G C 2 ' 1,Q: ([0, S 2 ) x M) x C 1 ([0,S 2 ),M s _ 1 ) in 
the limit i — > oo. Here, we use scaling invariant versions of the standard parabolic Holder 
norms, defined by ||w||co,a( M4 ) = ||«||c°(Af,s) + <5 Q Mc«(A/ a ) and more generally 

\\u\\* ca , b , HMs) = ^M^Wh^MsV M s = [0,S 2 )xM. 

k+2j<a 
j<b 

We remark that the second estimate of (A.l) immediately follows from Propostion 2.1 
and the Gronvall lemma, compare with (4.4). To estimate Wi = u.i — Ui-±, we observe 
that 

(A.2) d t w l - LiWi = fi 

for the elliptic linear operator 

LiW := A 9i w + A 9i (uj_i)(Vuj_i, Vw) + (dA gi (ui-i))(w)(Vui-i, Vit-s-i) 
and a right hand side that is bounded in C ' a (M$) for any S < Si by 

fi 2 \\fi\\c 'C(M s ) — C\\9i ~ 9 , i-l|lc ,1 ([0,5 2 ],Mf_ 1 ) + ^(IWIc' 2 > I ><»(M,s)) 

with a constant depending on a C 21,a bound on the previous iterate Ui-\ but not on Uj. 
We then apply the following scaling invariant version of parabolic Schauder estimates 

Proposition A.l. Let M be a closed manifold and let A > 0, A < oo be fixed. Then 
there exists a number C < oo such that the following holds true. Let L be any second 
order differential operator on Ms, 5 G (0, 1) any number, that is given in local coordinate 
charts as Lu = d Xi {a 13 d Xj u) + b l d Xi u + cu with 

a ij (x,t)^j > A |£| 2 for all £ G R m and (x,t) G M s 

\\°> ij \\h.°.°>(M s ) +Hb% a(Mt) +6 2 \\c\\* ca(Ms) < A. 
Then the solution u G C?> l > a {M s ) of d t w - Lw = f G C a {M s ) 7 w(0) = satisfies 

\\w\\* c *,i, HMs) <C5 2 \\f\\ 

C a (M 6 )- 

In terms of giving a proof of this result, we remark that standard Schauder estimates 
combined with a scaling argument give in a first step an estimate of the form 

(A.3) HIc=.i.«(m 5 ) < CS 2 \\f\\* ca{Ms) + C\\w\\ c o {Ms) 

for constants C independent of S. To retain this scaling invariance, we can then use the 
Ehrling-lemma and a further rescaling argument to estimate the second term of (A.3) by 

IH*)Hc°(M) < £ • IMIcWA/,) + C e Sup 6- l \\w(t)\\ L 3( Bs (x)) 

on every time slice. Finally considering the evolution of local energy quantities of the 
form J ip 2 \a 13 d Xi w{t)d Xj w{t) + S~ 2 w(t) 2 ~\ dx, ip a cut-off function supported on balls of 
radius 2S, gives that the last term in this estimate is bounded by a fixed (independent 
of S) multiple of <5 2 ||/||l=°, completing the proof of Propostion A.l. 

Turning back to the equation (A.2) satisfied by wt, this Schauder-estimate allows us to 
conclude that for any S < 5{ 

IMIcW.«(Af 5 ) ^ ct IK-i|IcW,«(m 5 ) + ^(IK-|lcW,a (Mj )) 2 . 
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Since u>j(0) = 0, the norm ||iw*||c73.i»<«(Afe) * s smau &t least for S small (a priori depending 
on i). We conclude that the first estimate of (A. 3) holds true, initially for 8 small 
and then, by a continuity argument, indeed for as long as the solution exists (provided 
= ^o(uo,ffo) was initially chosen small enough). But this very estimate prevents a 
blow-up before time <5j_i, so that Si = <5j_i = .. = So, completing the proof. 

A. 2. Proof of Lemma 2.8. We finally provide a possible proof of the fact that any 
metric in A4 S _ 1 can be written in the form f*g with / G T> s+l and g G M-\. 

Let s > 3 and let f2 C M a _\ be the subset of all metrics which can be written in the form 
g = f*g for a smooth metric g G -M_i and a diffcomorphism / of class H s+1 . 

We first prove that VI is an open subset of MLi using the slice theorem. Given any metric 
of the form g = ffigo, fo G V s+1 and go G M~\ we apply the slice-theorem 2.3 to the 
smooth metric go resulting in an A'llj-neighbourhood W of go, consisting only of metrics 
of the form f*gs, 9s an element of a slice S around go and thus in particular smooth. The 
pull-back W = fo*W is then an A^lj-neighbourhood of the original metric go G M s _i, 
containing only metrics of the form g = fo(f*9s) = (/ ° fo)*gs, Js £ S C M-i and 
/ o f G V s+1 . So indeed W <Zfl and tt is open. 

To see that fi is also closed, we use a result due to Ebin and Palais which says that 
the action of V s+1 on M. s _\ is proper, see e.g. Theorem 2.3.1 in [17]; in practice this 
means that if we are given a sequence of diffcomorphisms fi G T) s+1 and a convergent 
sequence of metrics gi — > g in Ai s _ 1 then knowing that f*gi — > g G M*_i converges (in 
H s topology) is enough to conclude that also (a subsequence of) the diffeomorphisms fi 
converges, fi—¥f in T> s+1 . 

Let now g G MLi be such that there are diffcomorphisms fi G T> s+1 and metrics gi G 
M-i with f*gi — > g (in M. s _i)- This convergence implies in particular that the length 
£(gi) = £(f*gi) of the shortest closed geodesic of {M,gi) is bounded away from zero. 
Thus the Mumford compactness theorem implies that after pulling back gi by a smooth 
family of diffcomorphisms /j, a subsequence of gi converges smoothly 

(fi)*9i = (/r 1 ° fiTiftm) ^ 9 G M-i- 

We conclude that the diffcomorphisms /r o fi converge to another diffeomorphism / G 
V s+1 and thus that g = (Z" 1 )*.^ G fi. 

References 

[1] W. Ding, J. Li, and Q. Liu. Evolution of minimal torus in ricmannian manifolds. Invent. Math., 
165:225-242, 2006. 

[2] W.-Y. Ding and G. Tian. Energy identity for a class of approximate harmonic maps from surfaces. 

Comm. Anal. Geom., 3:543-554, 1995. 
[3] J. Eells and H. J. Sampson. Harmonic mappings of riemannian manifolds. Amer. J. Math., 

86(1):109-160, 1964. 

[4] A. Freire. Uniqueness of the harmonic map flow from surfaces to general targets. Comm. Math. 

Helv., 70(l):310-338, 1995. 
[5] A. Freire. Uniqueness of the harmonic map flow in two dimensions. Calc. of Var., 3(1):95-105, 1995. 
[6] R,. D. Gulliver, R„ Osserman, and H. L.Roydon. A theorey of branched immersions of surfaces. 

Amer. J. Math., 95:750-812, 1973. 
[7] S. Kobayashi and K. Nomizu. Foundations of Differential Geometry. Interscience Publishers, 1963. 
[8] G. Licbcrman. Second order parabolic differential equations. World Scientific, 1996. 
[9] F. Lin and C. Wang. The analysis of harmonic maps and their heat flows. World Scientific, 2008. 
[10] M. Rupflin. An improved uniqueness result for the harmonic map flow in two dimensions. Calc. of 

Var., 33(3):329-341, 2008. 
[11] M. Rupflin and P. Topping. Flowing maps to minimal surfaces. arXiv:1205.6298vl. 



22 



MELANIE RUPFLIN 



[12] M. Rupflin, P. Topping, and M. Zhu. Decomposition of maps into minimal immersions, in prepara- 
tion. 

[13] J. Sacks and K. Uhlenbeck. Minimal immersions of closed riemann surfaces. Trans. Amer. Math. 
Soc, 271:639-652, 1982. 

[14] R. Schoen and S. T. Yau. Existence of incompressible minimal surfaces and the topology of three 
dimensional manifolds with non-negative scalar curvature. Annals of Math., 110:127-142, 1979. 

[15] M. Struwe. On the evolution of harmonic mappings of riemannian surfaces. Comment. Math. Helv., 
60:558-581, 1985. 

[16] P. Topping. Reverse bubbling and nonuniqucness in the harmonic map flow. Int. Math. Research 

Notices, 10:558-581, 2002. 
[17] Tromba. Teichmiiller theory in Riemannian geometry. Lectures in Mathematics ETH-Zurich. 

Birkhauser, 1992. 



Max-Planck-Intitut fur Gravitationsphysik, Am Muhlenberg 1, 
14476 Potsdam, Germany 



